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Abstract 

Polymer electrolytes typically exhibit diminished ionic conductivity due to the presence of cor- 
relation effects between the cations and anions. Microscopically, transient ionic aggregates, e.g. 
ion-pairs, ion-triplets or higher order ionic clusters, engender ionic correlations. Employing all- 
atom simulation of a model polymer electrolyte comprising of poly(ethylene oxide) and lithium 
iodide, the ionic correlations are explored through construction of elementary functions between 
pairs of the ionic species that qualitatively explains the spatio-temporal nature of these correla- 
tions. Furthermore, commencing from the exact Einstein-like equation describing the collective 
diffusivity of the ions in terms of the average diffusivity of the ions (i.e. the self terms) and the 
correlations from distinct pairs of ions, several phenomenological parameters are introduced to 
keep track of the simplification procedure that finally boils down to the recently proposed phe- 
nomenological model by Stolwijk-Obeidi (SO) [N. A. Stolwijk and S. Obeidi, Phys. Rev. Lett. 93, 
125901, 2004]. The approximation parameters, which can be retrieved from simulations, point to 
the necessity of additional information in order to fully describe the correlation effects apart from 
merely the fraction of ion-pairs which apparently accounts for the correlations originating from 
only the nearest neighbor structural correlations. These parameters are close to but not exactly 
unity as assumed in the SO model. Finally, as an application of the extended SO model one is able 
to estimate the dynamics of the free and non-free ions as well as their fractions from the knowledge 
of the single particle diffusivities and the collective diffusivity of the ions. 
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I. INTRODUCTION 



Polymer electrolytes have received increasing attention for electrolytic appli- 

cations in energy storing devices like batteries. Cations, ideally lithium ions (Li"^), supplied 
by a dissociated salt in a polymer medium provide for the charge conduction. It is known 
from experiments and computer simulations that the collective diffusivity (which is directly 
related to the ionic conductivity) in polymer electrolytes falls short of the average diffusivity 
of the charge carrying species (i.e. cations and anions). The reason is attributed to the 
presence of motional correlations between the unlike charges. 

Molecular dynamics simulation 5|, l6| of an amorphous polymer electrolyte offers the pos- 
sibility of exploring the nature of ionic correlations and their repercussions on the ionic 
conductivity. During the past decade realistic simulations on linear chains of poly(ethylene 



oxide) (PEO) with lithium salt have yielded insight into the ionic mobilities and conductiv- 



ities 
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,ll3|, |lj]. Typically, cations in polymer electrolytes can exist with or 
without the presence of anions in their vicinity. Cations which are coordinated solely by the 
ether oxygen (EO) atoms of the polymer host (i.e. PEO) are the so called free carriers and 
their proportion is important for the magnitude of ionic conductivity. However, a cation 
can also attract one, two or more anionic neighbors around itself forming the so called pairs, 
triplets or higher order clusters. The distribution of the oxygen and the anionic coordina- 
tion numbers is a consequence of the rich interplay between the interionic interactions and 
polymer entropy. Presence of ionic aggregates, specifically the neutral ones like the pairs, 
will not contribute to conductivity. 

The average diffusivity of the ions, Davg is related to the ionic conductivity, a, in an 
electroljd^e at temperature T by the equation 
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Here Ds and Xs are the diffusivity and the ionic fraction of species s, respectively; H the 
Haven ratio, Dcr 15(1 the collective diffusivity; and n is the total number of ions (both 



cations and anions), in a volume V . With H = 1 Eq. (ITal 
equation. D„ can be found from conductivity experiments 



expresses the Nernst- Einstein 
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be obtained from tracer diffusion measurements 
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17j or PFG-NMR 



19| whereas Dg can 
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2lj. Without 



ionic correlations the Haven ratio H is unity, yielding = Davg- If ionic correlations are 
present this relation needs to be augmented by the contributions from the cross correlation 
terms between the distinct ions, yielding H ^ 1. Correlated ionic motion can cause either 
H < 1 (as in inorganic ionic conductors) or H > 1 (as in polymer electrolytes). Note 
that in the latter case several cross correlations are present (cation-cation, anion-anion and 
cation-anion). 

Recently Stolwijk and Obeidi proposed a model (SO-model) that aims to explain the 
depression of ionic conductivity relative to the average ionic diffusivity in a representative 
electrolyte system of PEO and sodium iodide (Nal) IGj. This model is an attempt to deter- 
mine the mobilities of free ions and ion-pairs from the experimental information comprising 
of the tracer diffusivities of the ions and the ionic conductivity of the polymer electrolyte 
system. The key assumption of the model is that the ions exist either as free Na"*" and free 
1^ or in the form of contact ion-pairs (Na"^ ~ I )• Though the ion-pairs contribute to mass 
transport, they do not participate in charge transport due to charge neutrality, thereby, 
resulting in a diminished ionic conductivity. Their model predicts that, over a broad tem- 
perature range, the ion-pairs are orders of magnitude faster than either of the free ionic 



species 



16| . This led them to interpret that the ion-pairs are not coordinated with the ether 



oxygen atoms. In light of the insights from computer simulations where most of the cations 
(both free and those associated with anions) are structurally correlated and dynamically 
coupled to the ether oxygen atoms ^22^ the SO model can be extended to accommodate the 
scenario where the free ions and ion-pairs should display similar mobility. 

In this work we explore in detail the microscopic characteristics of the cross correlation 
contribution in a polymer melt via computer simulations for PEO with lithium iodide (Lil) 
as a salt. Our goal is fourfold. First, we derive an exact formal equation for D^r [viz. Eq. 
f l2b|) ] which fully contains the effect of dynamic correlations. Second, these correlations are 
interpreted in terms of simple physical pictures. For example it will be shown that ionic 
correlations exist beyond the nearest neighbor shell which implies that structural ion-pairs, 
i.e. those which are nearest neighbors, are not sufficient to grasp correlation effects. Third, 
we introduce a number of approximations such that the exact equation boils down to the 
phenomenological SO-model together with a few phenomenological parameters (all close 
to but not exactly one) by which we keep track of the approximations. For the present 
system, the values of these parameters can be retrieved from the simulations. This offers the 
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possibility of turning the SO-model more general and accurate. This attempt also brings 
out the necessity of additional information apart from merely the tracer diffusivities and 
ionic conductivity to quantify the microscopic dynamics. Fourth, we apply the extended 
SO-model to the experimental data, reported in [16]. It turns out that the data can be 
explained without invoking significant dynamical differences between the free ions and 
ion-pairs. The resulting fraction of associated ions from the model can be compared against 
the experimentally obtained information about the proportions of structural pairs via 
Raman scattering techniques. 



II. SIMULATION METHOD 

We have performed all-atom molecular dynamics simulation on a system of a poly- 
mer electrolyte comprising of chains of PEO and Lil as the salt in the canonical ensemble 
(NVT) using the GROMACS 231] package. The salt concentration in terms of the number 
of EO units and Li"*" ions in the system was EO:Li+=20:l. Conventional periodic boundary 
conditions were applied to get rid of surface effects. The simulation box contained 32 PEO 
chains. Each of the polymer chains had the chemical formula H(CH2— O— CH2)ArH where 

is the number of repeat units or chain length (Fig. [I]). In this work we have chosen 
N = 24 which is a reasonable chain length to study the dynamics of the ions/polymers 



within a production run time window o 



10 ns. The Nose-Hoover thermostat 



24 



5i 85 ns following an equlibration time of about 
251 ] algorithm was set to maintain an average tem- 



perature throughout the simulation run. The temperatures chosen were T = 425 K and T 
= 450 K. The relaxation time of the PEO chains in terms of the Rouse time, calculated 
by fitting the Rouse model 26| to the dynamics of the oxygen atoms were about 9.5 and 
4.3 ns at T = 425 K and T = 450 K, respectively. These relatively high temperatures were 
necessary to achieve equilibration and obtain reasonable estimates, in particular, for the 
collective diffusivity. Prior to the NVT runs the system density had been adjusted via an 
NPT (isobaric-isothermal ensemble with constant number of particles) run for durations of 
about 5 ns with an average pressure of 1 MPa. 

The starting configuration of the polymer melt for the NPT run was picked up from 
a thoroughly equilibrated system of neat PEO at a temperature of 500 K. The ions were 
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Figure 1: The atomistic model of poly(ethylene oxide) with the partial charges used are -0.0326 for 
carbon (cyan), -0.2792 for oxygen (red), 0.0861 for hydrogen (white) and -0.1187 for the terminal 
carbon atoms 

then randomly dispersed in it (i.e. positions of the ions were randomly selected) such that 
the minimum distance between a newly placed ion and the already existing particles (i.e. 
previously placed ions as well as atoms of the polymer) in the system was about 3 A. To 
obtain a stable co nfig uration, the system was then energy-minimised using the steepest 



descent technique [23j. The subsequent NPT runs produced converged average densities 
of 1082 gcm~^ and 1099 gcm~'^ at T = 450 K and T = 425 K respectively. The density 
fluctuations were of the order of 3.5 gcm~^. Following the NPT runs, the systems were 
further equilibrated at constant density and temperature (NVT) for more than 5 ns. In this 
paper we restrict ourselves to the T = 450 K data because of the better statistics. 

The force field for the description of intra- (bond lengths, bend-angle and dihedral 
potentials parameters) and inter-molecular (e.g. Buckingham potential parameters, partial 



charges) interactions of the 
ions are taken from Refs. 
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10 and the interactions between the PEO and the Li" 



271]. The interactions between the I ions and the PEO 
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and between the ions are essentially from Ref. |28| with slight modifications [29|]. The 
attractive part of the Buckingham potential between the I — Li, I — I and I — PEO have 
been slightly diminished. However, the charges (+le and — le for Li+ and I~, respectively) 
and the Coulomb interactions between the ions were unaltered and calculated using the 
particle-mesh Ewald technique within a distance cut-off of 10 A. This was done in order 
to generate a higher number of isolated ion-pairs and ion-triplets. However, the fraction 
of free and non-free ions in the system remained the same. Notwithstanding, with this 
modification the equilibrium polymer conformation of the sequence of dihedrals for the 
triad O— C— C— O was found to be in agreement with the published results of Ref. for 
their system of PEO/LiI with a composition of EO:Li = 15:1 at T = 450 K. The qualitative 
nature of the ionic correlations reported later in this article also turned out to be similar 
for both the original and the modified potential. In short, introducing the changes made 
to the potential only affected the state of ionic association without involving the average 
polymer conformations and dynamics. 



III. SELF DIFFUSION AND COLLECTIVE DIFFUSION 

The mean square displacement (MSD) of Li^ and I~ ions obtained from the MD simula- 
tion at T = 450 K are displayed in Fig. [2](a). The I~ ions are slightly more mobile (by a factor 
of ~ 1.6) compared to the Li"*" ions. The collective diplacement of the ions [see Eq. ([2])] which 
measures the mean square displacement of the center of charge of the system is also shown 
in Fig. [2](a). The collective diffusivity, Dcr is diminished by a factor of approximately 3 with 
respect to the average diffusivity of the ions, {0^ + Dj)/2. Fig.[2]^b) shows the experimental 
values of the tracer diffusivities of Na+ {D%^) and I (Df) in PEO/Nal (EO:Na=30:l and 
polymer chain lengths were of the order of 10^) from Ref. pL6|] as a function of temperature. 
The collective diffusivity for this system was reckoned from dc conductivity using Eq. ([1]). 
At T = 450 K factors of Df/D%^ ^ 1.7 and {D%^ + Df)/2D^ ^ 2.4 are observed for 
the experimental system similar to the simulated PEO/LiI system. The absolute values of 
the diffusivities of the ions in the simulated and the experimental system are of the same 
order of magnitude at the considered temperatures. Furthermore, the reduction of the ionic 
diffusivity [i.e. 0.5{Df + D^£-) — Da^] is close to the experimentally observed value. 
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Figure 2: (a) Mean square displacement of (MSD) Li^, and collective displacement from a MD 
simulation of PEO/LiI (EO:Li+ = 20:1) system at 450 K. The linear relationship between MSD 
and time is shown by a dot-dashed line, (b) Tracer diffusivities of Na^ i-^Na) ^ i-^f ) ^^'^ 
collective diffusivity (from dc conductivity experiment) in PEO/Nal with EO:Na+=30:l from Ref. 
[l^ . Note that in Ref. [3] the collective diffusivity D'^ was defined as D'^ = 2D„. 



From the simulation we identify a non-free ion as one which is associated with a counterion 
such that it is separated by a distance of less than 4.6 A that demarcates the first minimum 
of the Li — I radial distribution function, g^n [see Fig. [5] row 1, column 2]. Depending upon 
the number (one or two) of ions in the nearest neighbor shell of a Li+ ion one may speak 
of an ion-pair or an ion-triplet, respectively. The same cut-off radius is applied for defining 
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Table I: Proportion of free (n = 0) and non-free (n = 1: pairs, n = 2: triplets) ions in percentage. 



n 





1 


2 
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Li [I]n 


6.2 


75.9 


17.9 





I [Li]„ 


9.3 


69.8 


20.6 


0.3 




t[ns] 

Figure 3: Mean square displacement of free I~ (solid), free Li"*" (dashed), Li"^ with one I~ within 
its radius of 4.6 A (□) and Li"*" with two I~ ions within its radius of 4.6 A (dot-dash). 



a free or non-free I~. With this criterion the proportions of free and non-free ions present 
in the simulated system are displayed in Tab. [B 

We have computed the MSD of a subset of Li+ or ions which are either free or par- 
ticipating in an ion-pair or ion-triplet (see Fig. [3]). An ion was considered free during time 
t if it had no counterion within its nearest neighbor distance for at least 95% of the time 
t. Also shown are the MSD of those Li"*" ions which were paired with only one I~ (but 
not necessarily with the same I~) for 95% of the time. A similar definition was used for 
calculating the MSD of Li+ possessing two I~ neighbors (i.e. ion-triplets). This definition 
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limits the maximum values for t because of the finite life-time [see Ref. ll|]] of a specific 
local structure. The typical life-times of ion-pairs will be discussed in more detail further 
below. 

It is observed that the dynamics of the free Li"^ ions are almost indistinguishable from the 
paired Li"'' ions in dramatic contrast to the result of the data fitting, reported in 16|. The 
present result is expected since most of the Li+ ions in the simulated system are coordinated 
to three to six contiguous ether oxygen atoms of the PEO backbone and all cations are 
found to be complexed by EO atoms. Thus, the Li"*" dynamics mainly reflects the dynamics 
of the polymer chain. In particular, this strongly correlates the dynamics of free Li+ ions 
and paired Li+ ions. The dynamics of the subset of free I~ are faster than the free Li+ ions 
by a factor of 1.6. The MSD of Li+ ions in the triplets are marginally slower compared to 
the Li"*" in the pairs (Fig. [3]). The I~ in the ion-pairs are faster than the Li"*" in the ion-pairs 
for short times {t < 2 ns) and beyond this merge together. The difference at short times 
reflects the restriction of Li"*" motion due to the strong coupling to the polymer. 

Even though the typical life-time of a Li — 1 association only lasts for about 700 ps [Uj , 
there exist forward and backward jumps of the ions leading to an extension of the time 
duration that an ion is complexed to the same counterion. To see this clearly. Fig. H] shows 
the probability that a particular pair of Li"*" and I~ separated by less than 4 A exists at times 
t = to and t = to + t where to is some reference time. On an average about 6 ns is required 
for a Li^ ion to change its I~ neighborhood. Over this timescale a Li^ ion translates along a 
polymer chain by a distance of approximately 4 monomers. Thus, the Li"'" — 1^ pairs move 
collectively along a chain. This is seen from Fig. H] (inset) where the mean square variation 
of the ether oxygen index (An^(t)) experienced by a Li"^ while being still constrained to the 
same chain is plotted in time t (also see Ref. 22|). 
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Figure 4: Relative likelihood that a lithium ion is complexed to a specific iodide ion at times t = to 
and t = to + t. [Inset] Mean squared variation in the average ether oxygen index, (An^(t)), to 
which a Li*^ ion is coordinated while still being constrained to the same polymer chain in time t. 

IV. IONIC CORRELATIONS 



The collective diffusivity can be calculated from the microscopic dynamics of the ions via 



1 " 

= lim — V ZiZj{ARi{t)ARj{t)) (2a) 
t~^oo Qtn ^-^ 

^ n n 

= + E ^.^.(AR.(t)AR,(t))] (2b) 

i=l ijT^i 

= Davg - lim — — — (2c) 

^ t^oo Qt 

where ARj(t) = Rj(t) — Ri(0) is the displacement of a specific ion i carrying a charge of Zi 
during time t. Whereas the self terms [i.e. the first term of Eq. ( l2bl) ] express the diffusivities 
of all ions (i.e. Davg), the cross terms Scross{t) contain the correlation effects between distinct 
ions and give rise to a non-trivial Haven ratio. 

The second term of Eq. (I2cl) can be rewritten in terms of its different pair contributions. 
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Figure 5: Correlation functions: Each column displays a characteristic function of a specific pair 
of ions (e.g. Li — Li, Li — I, I — I) and each row displays the specific characteristic [e.g. radial dis- 
tribution function or rdf, gLAlir); directional correlation functions ALM{r,t), mobility correlation 
function BLM{r,t) and the correction term YiMif^t) in Eq. (jllh ]. The first row also shows the 
integrated rdf, | Jq dr Anr'^giMir) [dashed line]. 

Using ZLi+ = 1 and zj- = —1 one obtains 

—£-crossit) = £-LiLiif) — S-Lilif) + (3) 

where 

{L} {A/} 

^LM{t) = -'£Y.^^'^^^^)^'^^^t)). (4) 

Here, {L} denotes a set of cations or anions and the condition j ^ i'm the second summation 
of Eq. (jl]) implies that only distinct pairs of ions are considered. The negative sign in the 
middle term of the rhs of Eq. ([3]) results from the opposite signs of the charges on Li+ and 
ions. One expects in analogy to inorganic ion conductors a positive correlation of all 
ionic species, i.e. £^la/(^) > 0. However, since the correlation between Li"*" and I" ions will 
turn out to be particularly pronounced, for a polymer electrolyte, one gets Ecrossif) > 
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and H > 1 in contrast to inorganic ion conductors where one only has the contribution of 
cation-cation pairs. 

Our goal is to elucidate the origin of the correlations, contributing to ^^lm(^)- For this 
purpose we rewrite 

ARi{t,to)ARj{t,to) = ARi(t,to) ■ ARi(t,to) |ARi(t,to)| • |AR,(t,to)| (5) 

where ARj(i, io) = Ri(^ + ^o) — ^1(^0) denotes the displacement vector of ion i from 
its initial position at time t = to to its final position at t — to + t and ARj(i,io) = 
ARj(t,to)/|ARj(t, to)| refers to the unit vector of the displacement during time t of ion 
i starting from its initial location at t = to- Thus, a necessary condition for a non-zero 
^Luit) is the presence of directional correlations as expressed by the scalar product of the 
respective unit vectors. Qualitatively, one would expect that they mainly emerge from 
nearest-neighbor ions. Furthermore, the absolute value of £LM{t) might further increase, if 
nearby pairs are particularly mobile. 

To quantify these different contributions to SiMit) we introduce the following correlation 
functions between two distinct ions i and j 

aij{t; to) = ARi{t, to) ■ AR,(t, to) (6a) 
bij{t;to) = \ARi{t,to)\ ■ \ARj{t,to)\. (6b) 

The function aij{t;to) represents the directional correlation of two ions whereas bij{t]to) 
reflects their joint mobility during times t = to and t — tg + t. 

In a first step these functions will be averaged over pairs of ions separated by distance 
r. In general, the distance between two ions will change with time. A convenient way 
to take this into account is to use the initial and final inter-ionic separation distance with 
equal weight for the identification of r. With an additional average over to and using the 
abbrevation rij{t) = \Ri{t) — Rj{t)\ one can therefore define 

{aij{r, t)) = ^^{{Sr,rij{to) + Sr,ri,(to+t)}aij{t; to))to (7a) 

{bij{r, t)) = ^^{{Sr,rij{to) + Sr,rij{to+t)}bij{t; to))to (7b) 

where {■ ■ ■)to denotes the average over the time origins; Nij = {Sr,rij{to) + ^r,rij{to+t))to is 
the normalization constant [proportional to / drr^g{r)] and Sr,rij{t) is the Kronecker delta 
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function which bins the quantities aij{t; to) and %(t; to) at both the initial and final distance 
between ions i and j. Next, averaging the functions over all distinct pairs of ions, we define 



ALM{r,t) = {aij{r,t))i^{L};j€{M} 

BLM{r,t) = {bij{r,t))ie{L};je{M}- 



(8b) 



(8a) 



In the last step one can average ALM{r,t) over all distances. This yields 



^lm(^) = (2 - (5l,a/) / drnLM{r)ALM{r,t) 



(9) 



where we have defined 




(10) 



and p = n/V is the number density of the ions (cations plus anions). The prefactor 2 — 6l^m 
in Eq. ([9]) takes into account that the total number of unlike pairs (L ^ M) is twice the 
number of like (L = M) pair of ions . In Eq. (ITOl) the number density per species is expressed 
as p/2 and the radial distribution function (rdf) for the pair of species L — M is referred to 
as gLAiir). 

Likewise, one can define the quantity ELM{r,t) = {aij{r,t)bij{r,t))i^^L};je{M} using rela- 
tionships similar to Eq. ([7]) and ([8]). Our goal is to express ELhiit) (see Eq. H]) in terms 
of more elementary contributions. Based on the definitions, introduced so far, we can now 
write 



which implicitly defines the function liM(^;^)- If the directional correlation [via AiMi'^^i^)] 
and the mobility correlation [via BiM{''^-,iy\ of a pair of ions are statistically uncorrelated 
then one has i?LA/(r, t) = ALM{r,t)BLM{i",t), i.e. YLM{r,t) = 1. Possible correlations can 
be taken into account by slightly different values of YLM{r,t). Fig. [5] depicts the func- 
tions gLM^r), ALM{r,t), BLM{r,t) and YLM{r,t) for the different pairs of ionic species 
[(L, M) e {Li, I}]. The correlation functions, ALM{r,t) and BLM{r,t) are evaluated at 
specific times t = 0.1 ns, 1 ns and 4 ns. One can see that YiMif^jt) is larger than one. The 
resulting correlation between the directional and the mobility correlation, however, is rather 
weak. 



^LMit) = {2 - 6l,m) / drriLMir) ELM{r,t) 



(11) 



(2 - 6l,m) / drriLMir) YiMir, t)ALM{r, t)BLM{r, t) 
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Table II: Coordination numbers in the first, second and third neighborhoods using Eq. (jlOp . The 
bounds ri and r2 for each of the neighborhoods were determined from giAii^) in the Fig. [5j 



L - M 


"'LM 


(2) 


(3) 


Li - Li 


0.4 


3.8 


5.5 


Li - I 


1.1 


1.9 


3.7 


I - I 


0.3 


3.3 


5.2 



The high proportion of ion-pairs (Li — I) present in the system manifests as a pronounced 
nearest neighbor peak in the QLui^) at r = 2.9 A. In contrast, the rdfs between the like 
charges, gLiLi{f) and gii{r) are of similar strength in their respective first and second nearest 
neighborhood but exhibit weaker local density when compared to QLuir). To quantify this 
effect we introduce 

^Tm^^] drAnr^gLMir) (12) 

as the number of particles of species M in the k-th neighborhood of a particles of species 
L (Note that the coordination numbers of M around L and L around M, in the different 
shells, are the same, considering equimolar amounts of either ionic species). Here, r = ri 
and r = r2 are the minima of gLM{i^) corresponding to the k-th neighborhood. The values 
for n^^lj for A; = 1, 2, 3 are listed in Tab. [Tll 

In case that the system would have only free ions and well-defined Li — I pairs the 
quantity n^^^j would be identical to the fraction r„j of non-free ions. However, about 20 % 
of all ions form ion-triplets. This increases the value of n^^jj with respect to r„j. Numerical 
analysis shows r„/ ^ 0.92 (see Tab. [T]) which is, indeed somewhat smaller than n^ljj. The 
presence of these triplets also provides one important contribution to n^lin and nfj and thus 
to SLiLi{t) and Sii{t). Due to the smallness of n^lii^i and nfj one already anticipates that 
the contributions to Scross{t) from the like pairs is minor. 

The joint mobility functions, i.e. the average of the product of the length of the displace- 
ment vectors of two distinct ions during time t, BLM{f,t), are plotted in the third row of 
Fig. [5] for the species Li — Li, Li — I and I — 1. We see, that this function is essentially in- 
dependent of r for all times and is very similar for all the three types of pairs. In agreement 
with Fig. [3] this means that spatial proximity of two ions does not enhance its mobility. 

In contrast, ALM{r,t) decreases with r for all types of pairs. Qualitatively, this is due 
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Figure 6: Directional correlation contributions from the A;-tli neighborhood of B-type ions around 
A-type ions as a function of time, for the distinct pairs of ions of types {A, B) = {Li, I}. 

Contributions from first {k = 1; o), second {k = 2; □) and third {k = 3; a) neighborhoods are only 
displayed. 

to the screening of the interactions between the ions. Also, as time progresses ALAf(r, t) 
becomes smaller for fixed r. This is expected because each ionic aggregate has a finite lifetime 
beyond which it disintegrates. Each of the ions which were once part of an ionic cluster will, 
beyond a certain time, dislodge, migrate and finally participitate in the formation of another 
new cluster. We would like to mention that an ionic positional correlation function similar 
in spirit has been reported by Miiller-Plathe et al 0, S]- 

Among the general features of AlmIt, t) one can see that they decay with r in a step-like 
manner. The steps, representing constancy in the correlation strength, are short spanned 
and located at pair separation distances corresponding to the peak positions of the giKiij')- 
Thus, it may seem reasonable to characterize the average properties of the AiMi'^^ii) in the 
individual shells. For this purpose we define 

,{k) _ /;;c^r47rrVM(r)Az.M(r,t) 
j^,^ drA7rr^gLAi{r) 

where r = ri and r = r2 are the minima of 5'lm(^) corresponding to the k-th neighborhood 
(see Fig. [5]). Analogously one can define B^^l^{t) and E^^l^{t). In the following we analyze 
the first three shells ~ 1)2,3). Their time dependence is shown in Fig. [61 The 

functions are more or less constant at short times. However, these curves display a trend of 



decay with increasing time, as expected. Correlations for all pairs of ionic species show the 
order: A^llj{t) > A^^lj{t) > ^^^(t) where (L, M) e {Li, I}. Henceforth, we will often refer 
to a timescale t* = 0.5 ns. Note that t* is significantly smaller than the typical life-time 
of a Li - I pair. For t = t* we have A^^^jit") = 0.8, A^l^^^.{t*) = 0.53 and ^^^(r) = 0.49. 
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Figure 7: Parameters a, f3, Xi VLi-, fli (inset) measured as a function of time. 



Thus, the directional strength of motion of a pair of unlike charges (Li — I) is more potent 
in comparison to those of the like charges and is close to a perfect correlation A^^lj{t*) = 1. 
Furthermore, the shell-dependence is due to the screening of the higher shells of ions with 
respect to the central ion. 

After having discussed the individual contributions to ELMif) we now come back to Eq. 
(jl]) by rewriting it as 

SLMit) = (2 - 5l,m) l]nS,Efi(t). (14) 

k 

The correlations originating from the fc-th neighborhood is contained in the quantity 
E^lit). From the previous discussion we have seen that the dependence of E^l^l^{t) on the 
ionic pair L — M is mainly governed by A^l^\,j{t) which for Li — I is approximately 60 % 
larger than for the like pairs. Since, also, rS^Xj; is much larger for Li — I (see Tab. HT]) the 
by-far dominant contribution to £cross{t) stems from Smit). 
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V. STOLWIJK-OBEIDI MODEL 



In this section we begin with exact microscopic expressions and proceed through steps of 
simpUfications to construct an exact model similar in spirit to the phenomenological Stolwijk- 
Obeidi (SO) model. The specific procedure is determined by the input from simulations. The 
key idea of the SO model is to express the observed single-particles diffusivities, obtained via 
tracer diffusion, and the conductivity in terms of the properties of single ions and ion-pairs. 
In the SO model no triplets are taken into account. In case that these aggregates exist one 
should therefore better speak in terms of free ions or non-free ions, the latter possessing one 
or more counterion(s) in their coordination sphere, i.e. existing as both ion-pairs and ion- 
triplets in our case. In any event, this terminology gives rise to a more general applicability 
of that approach. 

First, we start to derive the single-particle dynamics. Here we have to take into account 
that during the life-time of a pair the MSD does not reach the linear regime. This requires 
a definition for the long time behavior of the transient ionic states. Specifically, we define 
the MSD associated with ions of type L (i.e. either Li+ or I^) existing in state S [i.e. either 
free (/) or non-free (nf)] in the following manner: 

/AP2 u\\ _ ^ie{L}i[^s{i,to),S + Ss{t,to+t),s]^^t{t,to))to 

\^-^L,s[t)) - 7T —J ^ ■ [^^) 

2^ie{L}\"<i,to),S + Os{i,to+t),s)to 

Here s{i,t) stands for the instantaneous state of the ion i, AR^(t,to) refers to the square of 
the displacement experienced by ion i between times t = to ^^nd t = to + t and the Kronecker 
delta function Ss{i,to),s picks the instantaneous state S of ion i at time to- For t <t* smaller 
than the typical timescale of the lifetime of nonfree ions this definition would give rise to the 
data in Fig. [31 For t 3> t*, (Ai?| j(t)) and (Ai?| „j(t)) start to merge because of averaging 
effects, i.e. free ions becoming non-free and vice versa. 
The fraction of free or non-free ions can be written as 

^^'^ = Tv^Ta Ta ^^^> 



Using this somewhat complex but generic definition one can show that Eq. ( ITSjl and ( fT6|l 
can be combined (see Appendix) to yield 

(Ai?i(t)) = rLj{^Rlj{t)) + rL,nf{^Rlnf{t))- (17) 



17 



Note that the motion of state specific MSD is only given for times smaller than the average 
lifetime of a pair of ions Li— I and thus in particular we choose t = t*. 
The average dynamics of a non-free ion is characterised by 

rLi,nf + ri^nf 

This definition can be used to relate the decrease in average ionic diffusivity to the correlated 
motion of the ions. Eq. ([3]) can be rewritten via the following steps at a chosen timescale 
t*: 

Scrossif) = Xi ^Lii{t*) (19a) 

= XiX2 n^'^E^llit') [cf. Eq. dH] (19b) 

= X1X2X3 n^'^ A^i},{e)B^ll{e) (19c) 

= X1X2X3X4 n^'^B'~ll{e) (19d) 

= X1X2X3X4X5 n^^\^Rlf{e)) [cf. Eq. m] (19e) 

= X1X2X3X4X5X6 rnf{^Rlf{t')). (19f) 



In the last step we have defined 



^nf = —■ (20) 



The parameters xi, X2, Xs^ X4) X5 Xe are the dimensionless correction factors required at 
every step of simplification in Eq. (fT9|) . Xi captures the cross correlations between the like 
charges. Because the missing contributions to Scross{t*) have opposite sign one naturally has 
Xi < 1. X2 accounts for ignoring correlations between Li+ and I~ ions beyond the nearest 
neighborhood (~ 4.6 A). One has X2 < 1 because of the presence of anticorrelations between 
the pair of ions beyond r > 16 A. The origin of the anticorrelations is not yet comprehended 
fully, however, physically this might arise from the tendency of the subensemble of particles 
to nullify its momentum due to a finite simulation box size. Xs is the correction due to 
the F-factor [see Eq. (fTTl) ]. which is given by the average value of YiMif^yt) in the nearest 
neighbor shell, refiecting correlations between the directional and the mobility correlations. 
X4 expresses the non-ideality of directional correlation in the nearest-neighbor shell, i.e. 
X4 ^ A^^^j{t*) < 1. The factor X5 is needed to approximate Bflj{t) by AR'^j^{t). If one 
compares the quantities AviArj and {Arf + Ar'j)/2 the latter is either greater than or equal 
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Table III: Parameters a, (3, 7, Xi^ VLi and rjj from simulation evaluated at time t = t* = 0.5. 



ait'') 


Pit*) 






X2it*) 


X3(t^) 


X4(t'^) 


X5(t^) 


X6(i*) 


VLiit*) 




0.959 


1.221 


0.966 


0.94 


0.88 


1.1 


0.86 


0.92 


1.18 


1.07 


0.82 



to the former. This rationahzes Xs < 1- Lastly in the final step the factor xe = n^^'^ /r. 



nf 



is 



required which as discussed above is larger than one. The individual values of x for t = t* 
are listed in Tab. IIIII 

Next we define the following parameters for later convenience: 



7 



Pit) 



ait) 



(21) 



ru,nf ■ (Ai?L,„/t))' {ARljit)) 

The time dependence of the parameters a, (3 and 7 is shown in Fig. [71 These are observed to 
be roughly stationary within the 1 ns time window which is less than the maximum life-time 
of an ion-pair. The values measured aX t = f* = 0.5 ns for the present system at T = 450 K 
can be found in Tab. Illli We have verified that choosing = 0.1 ns instead for defining the 
parameters would virtually yield similar magnitudes of results that are presented later. 

That the value of ait") is less than unity (0.96 and 0.91 for T = 450 K and 425 K 
respectively), can be understood from the slightly slower dynamics of the ion triplets than 
the free and the paired ions. It must be noted that a(t) ~ 1 is a consequence of the coupling 
between ion and segmental dynamics of PEO. On the other hand Pit*) ^ 1.2 means that 
the I~ ions in the pairs are slightly faster than their Li"*" counterparts at short timescales 
(see also Fig. [3]). It is interesting that the correction factor x(t*) is found to be less than 
unity (see Fig. [7]) unlike the assumption of xif) = 1 in the SO model. In essence, all the 
fixed parameters (Q;(t*), Pit*), 7(t*) and x(t*)) are close to unity due to the natural outcome 
of the ionic association and dynamics in the system and which is why the SO model could 
approximately define these characteristics. All these parameters are found to be nearly 
the same when going from the temperature of 450 K to 425 K (not shown). Even though 
Xi(t*) ■ ■ -Xeit*) differs slightly between the temperatures, the final product is found out to 
be similar due to compensation of the increase of one factor by the decrease of another factor 
and vice versa. 

The diffusion constant Dl can be written as Dl = lim^^oo(A_R|(r))/6r. Furthermore, 
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we define 

- }^ {ARlit^)) Qi (Ai?i(t*)) 

through which we are extrapolating the short time dynamics of the transient ionic state 
S of species L to long times using the time dependence of {R\{t)) as a reference. This 
extrapolation is necessary because {R\ si^)) '^'^y defined for times of the order of the 
lifetime of the ion in state S which is smaller than the timescale over which the ions become 
completely diffusive. As a specific advantage of the definition of D^ S) Eq. (fT7|l can now be 
rewritten as 

D'[ = TljDlJ + rL,nfDL,nf. (23) 

which is the form used in the SO model with L G {Li, /}. 

In the SO model all Xi ^i^re assumed to be one. According to the above discussion this 
corresponds to the case of ideal Li — I pairs and no directional correlations for other pairs 
(either pairs of like ions or pairs of Li"*" and I~ ions beyond the nearest neighbor shell). Mul- 
tiplying Eq. ([19|) with the factor £cross{t) /Qt£cross{t*) and using D^ross = \im.t^^£cross{t)/Qt 
one obtains from Eq. (119fl) and ( ITSl) 



2D cross = X rnf ( Zr^—r]LiD Li,nf + -^^i]iDi,nf) ■ (24) 
VI + 7 1 + 7 / 

Here, x is the product of the parameters Xi ' ' 'Xe- We have introduced the parameters rin 
and 7]! which are required for an exact relation between the diffusion coefficients in Eq. ( HM 
and are defined as 77^, = rjnit*) and rji = rjiit*) with 

This parameter i]L{t) (see inset of Fig. [7j) captures the difference in timescales over which 
{ARlj.^gg(t)) and (A_R|(t)) approach the linear regime, i.e. the respective ionic dynamics 
become diffusive. In the present case one has rj^i = 1.07 and rjj = 0.82 by choosing t* = 0.5 
ns at T = 450 K (see Fig. [7]). f]Li > Vi means that Li+ becomes diffusive later than I~ 
ions which is evident from Fig. [2]^a). As discussed in Ref. 22] the diffusive behavior of the 
Li"*" ions for long chains is achieved via jumps between different PEO chains. The residence 
time of the cation with one chain is found to be Tii ^ 100 ns [22]. Only for short chains (as 
in the present case) the timescale of Li"*" diffusion will be somewhat shorter because of the 
dominance of the center of mass diffusion of the polymer. In contrast, an I~ ion change its 
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cationic neighbor(s) quite frequently over a timescale of r/ ^ 7 ns and therefore is able to 
become diffusive faster than its cationic counterparts. 

In the following, we formally define a = ait*), {3 = P(t*), x = The set of three 

equations, [Eq. fl251) with L6{Li, 1} and Eq. fl24p ] is identical to the SO model if x = 1, 7 = 1, 
/3 = 1 and rjn = r]j = 1 i.e. if the picture of strict Li — 1 pairs would hold. However, these 
three equations contain four parameters (-Dlj,/, Du^nf, Djj, Vnf). Due to this ambiguity 
there are many possible solutions for a given set of parameters. The solution, presented in 
Ref. 161] involved additional assumptions about a similar temperature dependence of the 
parameters. The resulting solution was Dnj / Dn^nf ~ 1/74 for high temperatures. To 
circumvent this problem we have introduced the dimensionless parameter a{t) as in Eq. 
fl^T]) and determined its value (similar to 7(t), Xi{^)) from simulations. The reason for 
the observed magnitude of a{t) has been explained above. 

Using these parameters the relations of the generalized SO-model (see also Appendix II) 
can be solved: 

"■"^ XaiiLi{l+ h)D'[J Duress + a) ^ 

r)tr 

77/(1 + 7) 1 + 7 

The quantities D*-[^, Df and D cross are experimentally accessible. The parameters a, (3, 7, Xi 
rjLi and rji can be gathered from simulations, thereby permitting a solution for the unknowns 
r^/, Dii f and Djj from Eq. fl26l) . We neglect the possible temperature dependence of a, 
/3, 7, X arid TjL because it is difficult to find a priori arguments to estimate the temperature 
dependence of the four parameters. 



VI. APPLICATION OF STOLWIJK-OBEIDI MODEL 

First, we aply the SO model to our simulation data. We employ the parameters of 
Tab. [m] and D^^^ = 6.7 X 10^^ cmh-\ Df = 10.95 X 10"^ cm^s^^ and D^^ = 2.7 X 10^^ 
cm^s~^ (obtained through Fig. [2](a) from the simulation at T = 450 K) that gives D^ross = 
2.125 X 10^^ cm^s^^. These when plugged into Eq. (126!) yields the following: r„j = 0.915, 
D^ij = 6.966 X 10"^ cm^s-i and Dfj = 13.865 X 10"^ cm^s'^ and thus Djj/Duj = 2. 
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Table IV: Sets of parameters a, /3, 7, Xi VLi and rjj used for solving Eq. ([26 



SET 


a 


13 


7 


X 


VNa 


VI 


SET I 


0.959 


1.221 


0.966 


0.848 


1.065 


0.818 


SET II 


0.959 


1.221 


0.966 


0.848 


1.0 


1.0 


SET III 


0.959 


1.0 


1.0 


1.0 


1.0 


1.0 


SET IV 


74 


1.0 


1.0 


1.0 


1.0 


1.0 



D 



DLi.f VLx 



1.5. These, when compared 



From Eq. (1251) one should then obtain ^^^2^ 
to the actual values culled directly from the simulation, e.g., r„j = 0.92 (see Tab. [T]) and 
{AR'jj{t*))/{ARl^jit*)) = 1.6 (from Fig. [3]) are reproduced, in good agreement with the 
solutions of Eq. (!26l) . By letting rjn = rjj = 1 and the other parameters of Tab. Illll unchanged 
one obtains r^/ = 0.97 and Dij/Duj ~ 10 from Eq. ( !26|) . Clearly, the factors rjn and rji 
are important. 

Next, we apply the model on the same representative PEO/Nal polymer electrolyte with 
the concentration of EO:Na^ = 30:1 as Stolwijk and Obeidi did. Accordingly, we use the 

iracer diffusivities of Na+ (-D^^) and I~ (Df) and the collective diffusivity {D„) from Ref. 

IgI at different temperatures as input. Note the experimental data that have been used 
subsequently in this work is calculated back analytically from the fit results of Ref. IGj using 
the VTF temperature dependence of the single ion diffusivity prefactors and the Arrhenius 
temperature dependence of the pair formation constant (see details in Ref. {l6|). Four 
different sets for the parameters a, /?, 7, X; VNa and 777 are considered which we will refer, 
henceforth, as SET I to SET IV and tabulated in Tab. llVi This will allow delineation of 
the relative importance of the parameters. SET I has been retrieved from our simulation 
(same as in Tab. IIIII) and additionally, we assume this to be similar to the case of Na+ 
ion in PEO. SET IV corresponds to the original SO model of Ref. |16i]. Furthermore, we 
have superimposed a weak noise having a strength between ±5% to the experimental data 
[i.e D%^{T), D^[{T) and D(j{T)\ in order to evoke an impression about the stability of the 
solution (r„j, Dj^faj and Djj). The mean and standard deviation of the unknowns are 
calculated for each temperature from 100 independently noised realisations. We had also 
attempted mixing ±10% noise to the experimental data, however, the standard deviations 
of the solution exceeded mean values in the high temperature limit. Thus, the accuracy of 
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lO^/TfK"'] 10^/T[K^] 

Figure 8: Diffusivities of free I~ (solid line) and Na+ (dotted line) obtained from Eq. (I26p . Experi- 
mental data used is back-calculated from the final results in Ref. |ld| and then superimposed with 
a weak noise of ±5 %. The standard deviation is shown by the vertical bars. Parameters used for 
the different sets are listed under Tab. llVi 

the input experimental data is imperative in deriving the estimates of the solutions. 

Fig. El^a) shows the values of diffusivities of the free Na"*" (dotted line) and free I~ (solid 
line) obtained from the solution of Eq. (12^ with the parameters SET I. It can be observed 
that Djj is larger than D^aj over the temperature range shown. The curves monotonically 
decrease from the high to the low temperature end. Dij exceeds Djs/aj by about a factor 
of about 6—7 at the low temperature end to about 2—3.5 at the high temperature end. 
Particularly, at T = 450 K, Dij/D^aj = 3.6 which is about a factor of 1.8 larger than 
what is seen in the simulation (i.e. Dij/Duj ~ 2.0). However, this factor is expected to 
depend upon the chain length of the polymer host as shown in Ref. [2(|] in agreement with 
the current observations. By only supplanting the parameters rj^a = rji = 1.0 in SET I and 
letting the others unaltered (i.e. SET II) one observes a strong shift by a factor of ~ 3 in 
Dj J (dashed line) to larger values especially at the high temperatures and by about 20 % at 
the low temperatures. By choosing the parameters from SET III, Djj decreases less than 
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Figure 9: Fraction of non-free ions (rnf) versus inverse temperature (T) obtained from Eq. (j26p . 



16( 1 but superimposed with a weak noise 



Experimental data (back-calculated) used is from Ref. 
of ±5 %. The standard deviation is shown by the vertical bars. Parameters used for the different 
sets are listed under Tab. IIVI 

20 % at high temperatures and even less at low temperatures [see Fig. E^b)] when compared 
against SET II. Choosing the values a = 74 and /? = 7 = X = VNa = Vi = 1-0 SET IV) 
would have reproduced the results Dij/D^aj ~ 8 obtained by Stolwijk-Obeidi (not shown) 
at all temperatures. 

The dynamics of the free Na"*" ions exhibits insignificant alteration between the parameter 
sets I, II and III. This can be explained from Eq. ( l26|) (b) under the condition a ^ 1 which 
gives Dj^aj ~ D^N^ and therefore making Djyaj relatively independent of the parameters. 

In the foregoing, we have identified rin and 777 as important players in the estimation of 
Dj f while applying the extended SO model to our simulation data. Provided that a ~ 1, as 
is the case for the polymer electrolyte under investigation, a change in the other parameters 
(/?, 7 and x) are not as important in generating the large change in the diffusivity of the free 
I^ ions at the high temperatures. We, next, show that the reason for the large sensitivity 
of Dj f originates primarily from the small changes in 77/ whereas the influence of rj^a in 
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this case is rather weak. If a 1 then r/^„/ oc l/rj^a (from Eq. (l29l) and Eq. (l26l) (a)) and 
Dhnf oc rjNa/vi (see Appendix II) which implies that rj^nfDi^nf oc l/rji. From Eq. fl23l) . 
D/j = x-ri^;^'^^ 1 ^^"^ '"^.n/ large then the second term in the numerator of the rhs 
dominates and the denominator is also correspondingly small. Therefore, a small decrease 
in rii can cause a large increase in D/ j as observed in Fig. [^a) at high temperatures where 
the fraction of non-free ions is large (see below). In contrast, at lower temperatures the 
proportion of non-free I~ ions drops and thereby stabilizing Di j from small variations in 77/ . 

Fig. [9](a) shows the corresponding fraction of non-free ions for the sets I, II and III. 
Similar to the plots of diffusivities of the free ions, one finds monotonously decreasing and 
lightly bent behavior over the entire temperature range. Moving from SET I to SET II one 
observes an increase of r„j by about 6 % which can be explained from the proportionality: 
Tnf OC Xjrij^a- In coutrast, while changing parameters from SET I to SET III there is only 
a inconspicuous increase in the proportion of non-free ions by less than 2 %. This is due to 
the nullifying effect of the simultaneous variation of the parameters between the two sets of 
parameters and thus, appearing as similar fractions of non-free ions which is rather a mere 
coincidence. Choosing SET IV, however, results in a scanty fraction of non-free ions (see 

Fig.EKb)). 

Finally, the ion-pairing reaction constant, = t„j/(1 — t„j)^ [see Ref. 16|, l30|] for the 
equilibrium reaction Na+ + I^ ^ Na+I~ at the different temperatures is analyzed. Note 
that the above equation of kp is strictly applicable only if the non-free ions exist as ion-pairs. 
These are shown in Fig. for the four cases mentioned. For the original SO model (SET 
IV) the kp values are orders of magnitude less than those obtained from the other sets 
primarily due to the a parameter. A fit with kp = kpQexp{—AHp/kBT) where kpQ is the 
prefactor, AHp is the formation energy of the pair and fc^ is the Boltzmann constant yields 
the energetic term. This is listed in Tab. |V] along with the fraction of non-free ions at T 
= 450 K compiled for the different cases and sources. AHp calculated from the simulation 
using Vnf at T = 450 K is consistent with that obtained from the experiment by using 
the extended SO model with the parameters in SET I. Comparing the fraction of non-free 



ions between the original SO model 



16| (r„j = 0.08) against the extended SO model with 



SET I (r„j = 0.88) the discrepancy is evident at T = 450 K. In fact for the former a high 
ion-pair diffusivity is required to match the small fraction of ion-pairs in order to generate a 
reduction of the average ionic diffusivity by a factor of 3 to be compatible with the observed 
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Figure 10: Pair formation constant kp with parameters from the different sets Hsted under Tab. 

collective diffusivity. Note that with SET 11 (i.e. by neglecting the importance of rjL) one 
would have been tempted to erroneously conclude that the pair formation is not governed 
by a single AHp. As an additional information, the pair formation energy determined from 

the Raman spectroscopic work on poly (propylene oxide)PPO/NaCF3S03 with 0:Na = 30:1 

i 

by Kakihana et al. [SLf provides Aifp = 0.16 eV but with a high pair fraction. Of course, 
it would have made an interesting comparison if spectroscopic data had been available for 
the PEO/Nal system. 



I ' I ' I ' \ ' I ' ; 
'_ V --- SET II J 



— SET I 
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VII. SUMMARY 

In the first section ionic correlations in a polymer electrolyte (e.g. the archetypal 
PEO /Lil) stemming from ionic associations are elucidated by constructing elementary cor- 
relation functions. One of these functions included the time-evolution of the directional 
correlation of a pair of ions (both like-charged and unlike-charged) in dependence of the 
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Table V: Comparison of the energetic values obtained from fits of kp{T) from different somxes. 
The salt concentration in terms of the ratio of ether oxygen atoms to cations is included in the 
parenthesis. 



System 


Remark 


AHr, [eV] 


"•J 

(T = 450 K) 


PEO/LiI (20:1) 


Simulation 


0.31 


0.92 


PEO/Nal (30:1) 


Params, SET I 
(Tab. [ml) 


0.33 


0.88 


PEO/Nal (30:1) 


Ref. [30] 


0.19 


0.08 


PPO/NaCFaSOg (30:1) 


Ref. [31] 


0.16 


> 0.70 



separation distance between the ions. The correlations from the Li — I pairs are found to 
be the strongest in comparison to the Li — Li and I — I pairs of ions. Furthermore, the 
correlation effects are found to exist beyond the nearest neighborhood of counterions, albeit 
reduced due to the screening of interactions. In essence, the strong structural correlations 
together with the non-trivial values of the directional correlations between the cation-anion 
pairs lead to a reduction of the ionic conductivity of the polymer electrolyte compared to the 
average diffusivity of the ions. The joint mobility of a pair of ions calculated as the product 
of the magnitude of the displacement vectors of the ions showed that during a certain time 
this is almost independent of the structural correlations between the ions. One can thus 
infer that the dynamics of the different states of ions (i.e. free ions, ion-pairs, ion-triplets 
etc.) will be similar. Intuitively, due to the fact that most of the lithium ions are cou- 
pled to the ether oxygen atoms of the PEO one would expect some correlations between 
the dynamics of the free and the associated ions. This is, additionally, verified by comput- 
ing the mean square displacement of the free Li+, isolated Li+— I^ pairs and the isolated 
I^— Li+— 1~ triplets. This is in stark contrast to the conclusion from the Stolwijk-Obeidi 
model 16j which predicted the dynamics of ion-pairs to be orders of magnitude faster than 
the free ions. 

In the second section the SO model is extended by starting from the Einstein-like 
equation which describes the collective diffusivity of the ions from the single-particle 
dynamics of the ions (the self terms) and the correlations from the distinct pairs of ions. 
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Instead of describing the overall correlation effects only in terms of the dynamics and the 
fraction of isolated ion-pairs we generalised the approach by considering the dynamics 
and the fraction of non-free ions. Typically, the non-free ions account for the presence of 
ion-triplets or any higher order clusters along with the ion-pairs. Additionally, various 
approximation parameters, Xiy ^ire incorporated that keep track of the simplifications made 
to the correlated term of the Einstein-like equation. The simplifications can further be 
understood in terms of simple microscopic pictures. The parameters could be extracted 
from simulations and are found to be close to unity which is a consequence of the nature of 
ionic association and dynamics in polymer electrolytes. If each of these parameters except 
a (which appeared as a result in the SO model) are chosen to be exactly unity one recovers 
the SO model. The extended SO model along with the inference of dynamical similarity 
between the free and the non-free ions is applied to the same data as in Ref. 16| but with 



the additional parameters chosen from simulation. It is found that the fraction of non-free 
ions are high and consistent with the simulation results. The parameter 777- also turns out 
to be more important compared to the other parameters {rjNa, P, 7, x) specifically for 
yielding the estimates of the diffusivity of free I" ions. Typically, rjNa/rji 7^ 1 expresses the 
difference in timescales over which the cations and the anions become diffusive and reflects 
a characteristic of polymer electrolytes where the cations are predominantly complexed 
with atoms/groups located along the main backbone of the polymer for larger timescales 
in comparison to the anions which are coordinated to the same cationic neighbors for 
relatively short times. It must be, however, borne in mind that the results (i.e. fraction of 
non-free ions and the diffusivity of free ions) from the SO model depend on the accuracy of 
the input experimental data and possibly on the temperature dependence of one or more of 
the parameters (e.g. rjNa, rjj). Further work is required to map the values of the parameters 
as a function of interaction between ions and the resulting collective diffusivity of ions that 
would impart general applicability to quantifying ionic association. 



Appendix I 

Let n be the total number of ions in the system and so the number of cations and anions 
are n/2 respectively. Defining the quantity Z as follows 

n/2 

Z = X](^^(to),/ + ^sito+t)J + Ss(to),nf + Ssito+t),nf)to = ^ (27) 
i 
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one can show that Eq. (fT7|l can be derived from Eq. (fT5|) and Eq. (fT6l) through the following 
steps: 



z 

z z 

(28b) 



(28c) 



to 



s{to),nf + Ss{to+t),nf)to ^ '}2i{{^s(to),nf + ^s{to+t),nf}ARf {t))to (28d) 

Y.i{^s(tQ),nf + Ss{to+t),nf)to 

= Tu,f {ARljit)) + ru,nf (Ai?L,n/W)- (28e) 



Appendix II 

Here, we show how to rewrite Eq. fl2^ (for L = Li and L = I) and Eq. ([23]). The fraction 
of non-free Li+ ions can be determined as : r„j = 0.5(rij_„/ + r/ „j) = 0.5{rLinf + •yru^nf) = 
0.5rLj,„/(l + 7). Therefore, 

1 



and rj 



nf 



2r 



nf- 



7 



(29) 



+ 7 ^"^^ ■"1 + 7 
The corresponding proportion of free Li"*" and I~ ions are given by rnj = 1 — rn^nf and 
Tij = 1 — ri^nf- The solutions given in Eq. f l26l) can be obtained from the following equations: 



Du = rujDuj + rLi,nfaDLij = Dujiruj + arn^nf) 



D 



LiJ 



tr 



1 - (1 - a) 



1+7 



rijDij + aP—ri^nfDLij 



1 + 7 rjr 1 + 7 



2L). 



nj_ 

+ 7' 



(30a) 



(30b) 



(30c) 



2XTT^^L^(l+7/?)a^«,/ 
1 + 7 



where we have used Eq. f l2T]) and specifically in Eq. fl30al) we have used 
Di,nf = P^Dii^nf = CiP^Diij. Straightforward manipulation then gives rise to 
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Eq. ([26])(a-c). 
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